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Abstract 

We introduce an elegant method which allows the application 
of difTusing-wave spectroscopy (DWS) to nonergodic, solid- 
like samples. The method is based on the idea that light 
transmitted through a sandwich of two turbid cells can be 
considered ergodic even though only the second cell is er- 
godic. If absorption and/or leakage of light take place at the 
interface between the cells, we establish a so-called "multipli- 
cation rule" , which relates the intensity autocorrelation func- 
tion of light transmitted through the double-cell sandwich to 
the autocorrelation functions of individual cells by a simple 
multiplication. To test the proposed method, we perform a 
series of DWS experiments using colloidal gels as model non- 
ergodic media. Our experimental data are consistent with 
the theoretical predictions, allowing quantitative characteri- 
zation of nonergodic media and demonstrating the validity of 
the proposed technique. 



I. INTRODUCTION 

DifFusing-wave spectroscopy (DWS) [Q-^ is an exten- 
sion of "conventional" dynamic light scattering (DLS) 
to the multiple scattering regime. The basic idea 
of DWS is to use the autocorrelation function g2{T) = 
{I{t)I{t + T))j, I {I{t)fj, of the light intensity I{t) scat- 
tered by a turbid medium to study the dynamics of 
scatterers in the medium. Here (• • ■)rp denotes time av- 
eraging which can differ from the ensemble averaging 
{■■■)e- Introduced in 1987 |], the technique of DWS 
has rapidly evolved in recent years, and is currently ap- 
plied to study various types of turbid media, such as 



foams 
media. 



til 

It has 
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green bodies [|11|-[15[, emulsions 
granular ||2^ , |25| , and biological 
been demonstrated that DWS can be used to image 
macroscopic static and dynamic heterogeneities in turbid 
media |29|-35|. In 1995, Mason and Weitz have suggested 
that the motion of colloidal particles, characterized by 
DWS, can be directly related to the viscoelasticity of 
the surrounding medium (the corresponding experimen- 
tal technique is sometimes called "DWS-microrheology" ) 
p6| . Much attention has been paid to this approach 
since it could provide fast and non-invasive access to vis- 
coelastic properties of numerous materials, thereby open- 
ing a laxge new field of potential applications for DWS 

One of the main reasons for the remarkable success of 
DWS is the availability of a relatively simple and reli- 
able theoretical model which describes the experimental 
data in all practically important cases. The theoretical 
model of DWS is based on the diffusion approximation 
. It applies equally well to both statistically homoge- 
neous and heterogeneous media (media with inclusions, 
scatterer flows, etc.). Absorption of light, reflection of 
scattered waves at the sample boundaries, modulation of 
the source intensity, as well as various types of scatterer 
motion (Brownian |^,^ and sub-brownian motion, 
laminar [^-43| and turbulent flows, etc.) can be 
taken into account within the framework of the diffusion 
model. Recently, the theory of DWS has been extended 
to nonlinear random media [ ^J46| . 

An important condition for the applicability of the ex- 
isting diffusion theory to DWS experiments is the ergod- 
icity of the turbid medium under investigation. Indeed, 
ensemble-averaged quantities are commonly calculated 
theoretically, while it is the time averaging which is most 
easily obtained in experiments. Thus, (• • = {■ ■ 
is required for the experimental data to be described by 
the theory. If the light-scattering sample is nonergodic 
(say, the sample or some part of it is solid-like) addi- 
tional efforts, e.g. translational or rotational motion of 
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the sample during the measurement, are necessary in or- 



der to obtain (• ■ •)^ 



5[. Similar arguments hold 



for the role of nonergodicity in standard DLS experiments 

In the present paper we propose a novel, elegant and 
simple way of performing diffusing-wave spectroscopy of 
nonergodic media ||ll|,|5^. The method is based on the 
idea that light transmitted trough a sandwich of two tur- 
bid cells can be considered ergodic even if only the second 
cell is ergodic. We show that the resulting transmitted 
multiple-scattered intensity is ergodic despite the non- 
ergodicity of random medium in the first cell. Conse- 
quently, the usual DWS theory applies to the description 
of the intensity autocorrelation function 172 (t) . Moreover, 
the double-cell sample can be constructed in a way that 
the presence of the second, ergodic cell does not obscure 
the light-scattering signal of the nonergodic medium un- 
der investigation. In order to simplify the interpreta- 
tion of the experimental data, we introduce a so-called 
"multiplication rule" . Namely, we show that if a signifi- 
cant attenuation of diffuse waves occurs at the interface 
between the cells due to absorption and/or leakage of 
light in transverse directions, the ensemble-averaged /le/rf 
autocorrelation function gf^ (Li , L2, t) of the double- 
cell setup equals a product of autocorrelation functions 
g^\Li,T) and gj^'' (L2 , t) corresponding to the individ- 
ual cells [from here on, we use superscripts (1) and (2) 
to denote the autocorrelation functions corresponding to 
single- and double-cell samples, respectively; Li and L2 
are the thicknesses of the cells]. This allows a full DWS 
study of the turbid media in the first cell despite its non- 
ergodicity. Due to its simplicity and high statistical ac- 
curacy, our method is ideally suited for the study of tur- 
bid gels, colloidal glasses, and other nonergodic multiple- 
scattering media 15 5^ . The method also extends 
the applicability of the above-mentioned microrheology 
approach to an important class of solid-like com- 
plex materials where the motion of colloidal tracer par- 
ticles is highly constrained. 

The paper is organized as follows. In the next sec- 
tion, we provide a brief review of known experiments in 
the field of DLS by nonergodic media. We describe the 
methods used to overcome the problem of nonergodic- 
ity, and discuss the applicability of similar methods in 
DWS experiments. In the second part of the section, we 
introduce an original, novel method to deal with noner- 
godicity in DWS and reveal the conditions under which 
the method can be applied. Section III is devoted to 



the theoretical model of DWS in a double-layer medium. 
Starting from the diffusion equation for the field autocor- 
relation function, we obtain an expression for the auto- 
correlation function of light transmitted through a sand- 
wich of two turbid cells, separated by a nonscattering 
but perhaps absorbing wall. Section IV introduces the 
so-called multiplication rule, and reveals the conditions 
of its applicability. In Sec. ^ we describe our experi- 
mental setup, and Sec. VI is devoted to the discussion of 



the main results of our experiments. First, we have per- 
formed model experiments to test the validity of our the- 
oretical model, and to ensure that the experimental setup 
is adequately described by the theory. In these experi- 
ments, no attenuation occurs on passage of light through 
the wall separating the light-scattering cells, and the mul- 
tiplication rule does not apply. Next, we show how our 
method can be applied to an important case of consider- 
able leakage of light at the inter-cell wall. This situation 
is the most interesting for practical applications, since 
the multiplication rule implies that the autocorrelation 
function g'-^\Li,T) corresponding to the first layer can 
be obtained simply by dividing gl ' {Li, L2,t) measured 

for the two-cell setup by g'^\L2,T) measured for the sec- 
ond cell taken alone. In addition, our second sample is 
highly asymmetric: the first layer is of significant opti- 
cal thickness and has a relatively short correlation decay 
time Ti , while the optical thickness of the second one is 
moderate and its correlation decay time T2 S> ti. This 
prevents the light-scattering signal of the first, noner- 
godi c ce ll from being obscured by the second cell. In 
Sec. VII we discuss different aspects of application and 
optimization of the two-cell technique in practice. We 



summarize our main results in Sec. VIII 



II. 



DYNAMIC LIGHT SCATTERING BY 
NONERGODIC MEDIA 



In nonergodic random media, the scatterers are local- 
ized near fixed average positions, probing only a small 
fraction of their possible spatial configurations by ther- 
mal motion. As a consequence, the measured time- 
averaged quantities (such as the scattered intensity or 
its autocorrelation function) differ from the ensemble- 
averaged ones. Experimentally, one finds that a series of 
(time-averaged) measurements on a given sample yields 
a set of different results, each being of limited use for the 
characterization of the medium (see Fig. |l|). 



A. Concepts to deal with nonergodicity in dynamic 
light scattering 

For diluted, nonergodic samples several methods have 
been applied to properly average the data obtained in 
DLS experiments [[47|-p3[ . A good comparative review of 
some of these methods can be found in Refs. Q and gl) • 
The most direct method of performing ensemble averag- 
ing of scattered light for nonergodic samples is based on 
the idea of collecting light scattered by different parts of 
the sample, thus performing the "real" ensemble averag- 
ing pO[ . Experimentally, the sample is slowly moved or 
rotated, while the autocorrelation function of the scat- 
tered intensity 52(1') is collected. 
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FIG. 1. Typical light scattering signal of a non-ergodic sys- 
tem. Repeated (time averaged) measurements of g2{T) — l lead 
to non-reproducible results, each being of limited use for the 
characterization of the system. Over the measurement time 
(typically a few minutes), the time and ensemble averages of 
the fluctuations of the scattered light intensity are not the 
same. 

Obviously, this leads to an additional decay of 92(^)1 
which becomes increasingly rapid with increase of the 
translation velocity or the rotation frequency |50[| . The 
rotation/translation method can be extended to concen- 
trated, turbid suspensions without any particular difficul- 
ties, and it has been actually em ployed in recent experi- 
ments on diffusing- wave imaging ||3j,^. It has, however, 
an important disadvantage of experimental complexity 
(translation or rotation of the sample is required). Be- 
sides, it is not well suited for fragile turbid systems (such 
as gels), since moving or rotating the sample can lead to 
its disturbance or even breakdown. 

Pusey and Van Megen have proposed a method 
of obtaining the ensemble-averaged intensity autocorre- 
lation function from a combination of static and dynamic 
experiments. The idea is to measure the autocorrelation 
function of scattered intensity for a single sample orien- 
tation, and to subsequently rotate/translate the sample 
rapidly in order to obtain the average intensity for a given 
scattering angle [^7|- ^^ . Recently, an extension of the 
above method to DWS experiments has been developed 

ii- 

In general, all these approaches rely on a one dimen- 
sional motion of the sample, i.e. scanning the speckle 
pattern, which is not very efficient and therefore time 
consuming. This makes it intrinsically difficult to ap- 
ply these methods to systems which evolve in time, e.g. 
gelling systems. 

A different method of dealing with nonergodicity in 
DLS has been proposed in Refs. |Q and Q (see also 
Ref. iQ). The authors use a CCD camera to record 
the temporal evolution of many speckle spots simultane- 
ously, which allows them to perform correctly both time 
and ensemble averages (the latter is the average over a 
large number of speckle spots). The method appears 
to be very efficient in the small-angle single-scattering 



regime, while its application to multiple-scattering sys- 
tems is complicated by a (generally) insufficient time 
resolution of available CCD cameras (in DWS, fast and 
low- intensity speckle fluctuations are usually monitored). 
However, the method of Refs. and can be ap- 
plied to study extremely slow dynamics in turbid systems 
where the above mentioned constraints do not apply [ p9| . 

B. Two-cell technique 

To overcome the problem of nonergodicity, we pre- 
pare a sandwich consisting of two independent glass cells. 
The first cell contains a solid-like nonergodic medium un- 
der study, while the second cell is filled with an ergodic 
medium [see Fig. ^(a)]. By adjusting the concentration 
of scatterers in the second cell, the viscosity of the liquid 
where the scatterers are suspended, and the thickness of 
the cell, it is possible to shift the "forced" decay of the 
autocorrelation function due to the second cell to long 
correlation times t. Then g\ (Li,L2,t) will exhibit a 
short-time decay due to the motion of scatterers in the 
first cell, at intermediate r it will saturate at a plateau be- 
cause of the nonergodic nature of the medium in the first 

(2) 

cell, and finally, at long correlation times, g\ (Li,L2,t) 
will decrease to zero due to the motion of scatterers in 
the second cell. Physically, the slow motion of scatter- 
ers in the second cell gently shakes and randomizes the 
speckle pattern of the nonergodic medium. In this way, 
ensemble averaging of light scattered by the first cell is 
accomplished, similarly to the case when the sample is 
moved. The time scale of this averaging can be well con- 
trolled by the scatterer motion inside the second cell, 
and by the thickness of the cell. The only constraint be- 
ing a high enough optical thickness of the second cell, 
ensuring that the nonergodic intensity fluctuations pro- 
duced by the first cell are averaged out upon the trans- 
mission through the second one, and that no light can 
pass through the latter without being scattered. The de- 
scribed setup offers an advantage of high statistical ac- 
curacy due to a two dimensional averaging scheme. Fur- 
thermore, no mechanical disturbance is applied to the 
sample, which is of particular importance when working 
with fragile systems. 

Let us briefly discuss different aspects of time and en- 
semble averages in the double-cell geometry. For sim- 
plicity, we choose the intensity of light (and not its auto- 
correlation function) as an object of averaging, but the 
same arguments apply to the autocorrelation function 
as well. First, we consider the transmission of a plane 
wave through a layer of totally rigid random medium 
(thickness Li, photon transport mean free path Z*)- We 
flnd that the ensemble-averaged intensity of multiple- 
scattered light (/(R, z))^ at depth z and transverse po- 
sition R = {x, y} is independent of R. Meanwhile, 
(/(R, z))rp = /(R, z) is a random function of R, a "fin- 
gerprint" of a given scatterer distribution in the sample. 
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Nonergodic medium : 
scatterer motion is restricted 



Nonscattering wall 



Ergodic medium: 
scatterers move freely 
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Incident intensity profile 
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between tiie layers 
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Time-averaged transmitted intensity ^ 



FIG. 2. Schematic illustration of the idea of the two-cell 
technique, (a) A cell (thickness L2) containing an ergodic 
turbid medium is added just after the cell containing the non- 
ergodic medium under investigation (thickness Li). The cells 
are separated by a glass wall (thickness A). The overall thick- 
ness of the sample is L = Li -f A -I- L2. (b) Time-averaged 
intensity profiles are shown schematically for the incident 
light, the multiple-scattered light at the inter-cell interface 
z = Li, and for the transmitted light. Between the cells, the 
time-averaged intensity exhibits speckles due to the noner- 
godic nature of the medium in the first cell. The speckles are 
washed out upon the transmission through a sufficiently thick 
second cell. 



port is valid in both layers. This allows us to estimate 
the typical speckle size at the z = Li plane as d ~ A 
(see, e.g., Ref. [^), where the wavelength of Hght A 
is assumed to be roughly the same in both layers and 
much smaller than the scattering lengths in the layers. 
At the same time, as long as L2 3> ^J, the average 
intensity of light transmitted through the second layer 
is not sensitive to those details of the "source distribu- 
tion" /(R, Li,t) which are finer than L2. Consequently, 
the speckle pattern /(R, Li,t) will be completely washed 
out upon the transmission through the second layer, and 
(/(R, L, t))rp will be insensitive to the exact spatial con- 
figuration of scatterers in the first (rigid) layer, being 
equal to (/(R, L, t))^. The scattering system is then er- 
godic: the time and ensemble averages are equivalent. 
Similar arguments apply in the case when the first layer 
is not completely rigid, but exhibits some nonergodic dy- 
namics. The general reason for this is the following. A 
given speckle spot of typical size d <C ^2 the z — Li 
plane can be considered as a point source of light, which 
produces a photon cloud spreading linearly upon the dif- 
fusion through the second layer of turbid medium. Hence, 
the speckle spot at z = Li gives contributions to the av- 
erage intensity within a region of typical transverse size 
L2 at the z — L plane. For this reason, the intensity of 
light at a given point R of the z — L plane is formed as a 
sum of contributions of '--^ (^2/^)^ S> 1 speckle spots 
located at z = ii [see Fig. ||(b)]. In this way, the aver- 
aging over a large number of independent speckle spots 
is performed. In some sense, such an averaging is equiv- 
alent to that obtained by translation (or rotation) of a 
single nonergodic layer. It is worthnoting, however, that 
in order to average over the same number N of speckle 
spots, the sample has to be translated by a distance of 
order Ll/d'^ L2- 



III. THEORY OF DWS IN A DOUBLE-LAYER 
MEDIUM 



Consequently, (/(R, z))^ / (/(R, z;))^, and the consid- 
ered light-scattering system is nonergodic. 

Let us now add a second turbid layer (thickness ^2, 
photon transport mean free path ^2) just after the first 
one [we neglect for a moment the thickness A of a glass 
wall separating the cells, see Fig. ||(a)]. We assume that 
the scatterers in the second layer are moving, and that all 
possible spatial configurations of scatterers are explored 
in course of this motion. The light intensity /(R, z, t) 
is then time-dependent, and the time- averaged intensity 
of light transmitted through the double-layer medium, 
(/(R, L, t))j,, can be considered as a result of transmis- 
sion of the speckle pattern (/(R, Li, t)) j,, being specific 
for the given configuration of scatterers inside the first 
layer, through the second layer. Now we assume that 
the optical thicknesses of both layers are large enough 
(L„ ^ /*), so that the diffusion picture of wave trans- 



Adding a second cell modifies the intensity autocorre- 
lation function of transmitted light, which now exhibits 
an additional decay. In order to interpret properly the 
experimental data, we need a theoretical model describ- 
ing the autocorrelation function in a double-layer turbid 
medium. In this section, we show that such a model can 
be readily constructed, provided that the transport of 
light is diffusive in both layers. 

Consider a slab situated between the planes z = 0, 
z = L = Li+A + i2, and consisting of two lay- 
ers of turbid media (thicknesses Li and L2) separated 
by a nonscattering but perhaps absorbing wall (thick- 
ness A, absorption coefficient Ma), as depicted in Fig. 
||(a). Li and L2, are assumed to be much larger than 
the photon transport mean free paths inside the layers, 
l\ and respectively. If the slab is illuminated by a 
plane monochromatic wave, the field correlation func- 
tion G'i(r, r) = (-E(r, t)E*{r, t + t)) ^ inside each turbid 
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layer can be found as a solution of the diffusion equation 
1,1: 



[V^-a^ir)] Gi(r,r) 



35(r) 

I* '' 



(1) 



where Q!^(t) describes "attenuation" of correlation due 
to scatterer motion, and 5(r) — sq5[z — zo) is the source 
term (zq ~ Z*). The explicit form of the "attenuation" 
coefficient cP'{t) in Eq. (|l|) depends on the type of scat- 
terer motion (e.g., Brownian motion, flow, etc.). One can 
show that (r) is determined by the absorption length la 
and the mean square scatterer displacement {|Ar(r)^)^ 

"lilJer 
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TJ* 



1*2 



(2) 



where k = 27r/A is the wave number of light in the 
medium. For Brownian motion |^,^, (Ar(r)^)^ — 
(Ar(T)2)^ = QDbt, and a^ir) = i/ilj*) + Qt/^t^I*'') 
with To = {k^DB)~^, while for a directed flow one finds 
(Ar(T)2)^ ^ oc r^, and a2(r) cx In the case of 
a gel, which is an example of nonergodic medium consid- 
ered in this paper, we assume that the particles undergo 
a sort of arrested subdiffusive motion |fl]l,p8^,p6| : 



(Ar(r)2) 



E.T 



I — exp 



(3) 



where S denotes the maximum r.m.s. particle displace- 
ment, Tc is a characteristic time required for this displace- 
ment to occur, and p is a free parameter. We expect fll]] 
p ~ 0.7 ± 0.1. Eqs. (|) and (|) give 



{k5f 
1*2 



1 — exp 



(4) 



The time autocorrelation function for a sample con- 
sisting of two turbid layers separated by a nonscattering 
wall can be found by solving Eq. (|l|) inside each layer, 
and then applying the boundary conditions at the sur- 
faces of the medium (z — and z — L) and at the in- 
terfaces between the turbid layers and the nonscattering 
wall. Boundary conditions at z = and z = L are of a 
well-known form 



[Gi(r,T) - zi (e, • VGi(r,T))] U^o = 0, (5) 
[Gi(r,T)+Z2(e. ■VGi(r,T))]|,=L -0, (6) 

where is a unit vector parallel to the z-axis, 2„ = 
(2/3)Z*, and the refractive index mismatch between the 
scattering and transparent media is neglected. In gen- 
eral, a rigorous theory would require using the extrapo- 
lation lengths Zn calculated with account for the refrac- 
tive index mismatch at the sample surfaces ]60|-|63[|, or 
deduced from the angular distribution of diffusely trans- 
mitted light [ p4| . For our purposes, however, it is suf- 
ficient to know that z„ ~ since the actual values of 



Zn are of no importance in the limit of ctn/* 1 and 
Z*/L„ ^ 1 that we apply in the following. 

Boundary conditions at the interfaces z — Li and 
z = Li -|- A between turbid and nonscattering media are 
found by applying the condition of flux conservation, as 
discussed by Ripoll et. al. Neglecting the refractive 
index mismatch between the turbid and nonscattering 
media, we find [[7l| : 

VGi(r,r))]U.L, 

+ Z2(e, • VGi(r,T))]|,=i,+A , (7) 
VGi(r,r))]U=i,+A 
-zi(e,. VGi(r,T))]|,=i, , (8) 

with /(r) describing the losses of energy at the inter- 
layer wall due to absorption and/or leakage of light in 
transverse directions [0 < /(r) < 1 and /(r) = 1 in the 
absence of both absorption and leakage] . For an infinitely 
wide slab (no leakage) we have |7^ 

/ = exp(-AfaA)(l - Ma A) + (M,A)2r(0, M„A), (9) 

where r(. . .) is the incomplete gamma function. For real- 
istic samples of finite width, leakage of light in transverse 
directions may be considerable. In the absence of absorp- 
tion {Ma ~ 0), we find for a cylindrical sample of radius 
R ^ l!^ after averaging over the sample crossection: 



[Gi(r, r) + zi (e^ • 
= /(r) [Gi(r,r) 

[Gi(r, r) - Z2 (e^ • 
= /(r) [Gi(r,r) 



4+ = 



1/2 



(10) 



Averaging over the sample crossection makes / indepen- 
dent of r, which largely simplifies the further analysis 
without affecting the final result qualitatively. 

In what follows, we assume that /(r) in Eqs. (|^) and 
(^ is given either by Eq. (||) or by Eq. (^, being in- 
dependent of r in both cases. The case when both ab- 
sorption and leakage are present can also be analyzed 
without any particular difficulties, but this leads to cum- 
bersome formulas without introducing any new qualita- 
tive features. Both Eqs. (^) and ( p^ ) exhibit a monotonic 
decrease from 1 (M^A = or A/R = 0) to (Af^A > 1 
or A/R ^1). If / = 1 (no absorption, no leakage), 
Eqs. (0) and (||) reduce to the well-known boundary 
conditions between two turbid media in a direct contact 
p9| , ^|32| . If / = (strong absorption and/or leakage), 
Eqs. (|7|) and (^) decouple and the problem is reduced to 
two single-layer problems. 

Equations (|l|) and (|1)-(H) enables us to calculate the 
autocorrelation function of light transmitted through the 
double-layer system, Gi{L, t)/Gi{L, 0), which in the ab- 
sence of absorption {la — > oo) inside the turbid layers 
reads 



5f^(Li,L2,T) 



aiLi 



02^2 



sinh(Q;iii) sinh(Q!2-/j2) 



X (l-/2 + 2zi/Li + 2z2/i2) 



1-/2 + 



2q:iZi 



2q;2Z2 



tanh(Q;iii) tanh(Q;2L2) 



(11) 
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where we assumed q;„Z* ^ 1 and 1^1 1. 

In our experiments, the time-averaged autocorrela- 
tion function of the scattered intensity (72(1", t) = 



{I{v,t)I{v,t + r))^/{I{v,t% 



measured. As the 



double-layer system is shown to be ergodic (see Sec. ||), 
(72 is related to 171 by the Siegert relation 



52 (r, 



l + /3|5i(r,r)| 



(12) 



where /3 ~ 1 is a constant determined by the experi- 
mental setup 1^,^. Let us assume that the first layer is 
filled with a gel, while the second one contains a sus- 
pension of Brownian particles with a correlation decay 
time T2 = To{l2/L2)'^- We also assume that the char- 
acteristic correlation decay time of the first layer ti is 
much smaller than T2. Then, it can be shown from Eqs. 
( pi] ) and ( p2[ ) that when the r.m.s. particle displacement 

(Ar(T)^)^ ^ approaches 6 for r ti [Eq. (||)], the auto- 

correlation function g\ {Li, L2,t) reaches a plateau for 
Ti ^ T <C T2. The plateau height can be found from 
Eq. (|rT|). As T approaches T2, g^j^\Li, L2,t) continues 
to decrease. 



IV. MULTIPLICATION RULE 

The autocorrelation function of light transmitted 
through a single layer of turbid medium can be found 
PH from Eq. (|l|) with boundary conditions (||) and (||), 
and since for al* <^ 1 and L I* the solution reads 



aL 



sinh(aZ/) 



(13) 



We can rewrite Eq. ( |11| ) as 

g^;'\L^,L2,r)=g['\L^,T)xg['\L2,r)xF, (14) 
where 



F={l-f + 
1-/2 + 



2zi/Li + 2Z2/L2) 
2aiZi 



2a2Z2 



tanh(Q!iLi) tanh(Q;2i2) 



(15) 



Hence, the autocorrelation function of light transmitted 
through the double-layer sample is given by a product of 
autocorrelation functions of individual layers, times some 
function F which describes the coupling between the lay- 
ers. By adjusting the parameters of the experimental 
setup, ~ 1 can be achieved, and then the following 
multiplication rule will hold: 



gP(Li,L2,r)^.g(^^(ii,r)xgi^^(L2,r). 



(16) 



To reveal the conditions of validity of the multiplication 
rule (|l|), we consider separately the cases of low (/ ~ 1) 
and considerable (/ < 1) losses of light at the inter-layer 
interface. 



If / ~ 1, Eq. dH) becomes 



F : 



£1 
'1 
Li 



L2 

I* 
'2 

02^2 



tanh(Q;2i2) I2 tanh(aiii) 



(17) 



which reduces to unity only if aiLi <C 1 and a2^2 <C 1, 
or Li/ll :s> L2/I2 and a2i2 ^ 1, or Li/l'l ^ ^2/^2 and 
aiLi <C 1. In other words, for the multiplication rule 
( pj| ) to hold the optical thickness of e.g. the first layer 
should be much greater than that of the second one, while 
g[^\L2,T) ~ 1 for the latter. 

In the presence of losses at the inter-layer interface, 
/ < 1 and Eq. (^ yields F ~ 1 if i„/Z; > 1/(1-/^) and 
Q!„Z*/ tanh(Q;„L„) ^ 1— Z^. To give an example, suppose 
that 50% of wave energy is lost on crossing the inter- 
layer interface once (/ ~ 0.5). This gives Ln/ln ^ 4/3 
and a„Z*/ tanh(a„L„) 3/4. The latter conditions are 
commonly satisfied in typical DWS experiments, making 
the application of the multiplication rule rather practical. 

It is worthwhile to note that there exists a different, 
less rigorous, but more transparent way of establish- 
ing the multiplication rule. Adopting the path-integral 
picture of light propagation through the double-layer 
medium , we can write the field autocorrelation func- 
tion of transmitted light as 

/"OO /-OO 

gf ■'(Li,L2,t) = / dsi / ds2P2(si,S2) 
Jo Ja 



X exp 



(18) 



where P2(si,S2) is the relative weight of paths consist- 
ing of segments of overall lengths si and S2 inside the 
first and the second layers, respectively. The variances 
of phase differences ^A(p2(r))_^ describe the dephasing 
of light inside the layers. For a single layer of thickness 
L one has 



g^^\L,T) = / dsFi(s)exp 



(19) 



Obviously, if the variables si and S2 can be considered 
independent, ^2(51,52) = Pi{si)Pi{s2), and Eq. ( p^ ) re- 
duces to a product of two terms, each being of the form 
(p9|). This transparent physical picture underlying the 
multiplication rule is particularly helpful if one seeks to 
understand the role of losses (absorption and/or leakage 
of light) at the inter-layer interface. Increasing the losses 
makes the partial lengths of light paths si and S2 in the 
layers more and more independent, since the losses reduce 
the probability for a typical path to cross the interface 
more than once. Consequently, most of the paths will 
consist of two independent segments, one inside the first 
layer (length si), and the other one — inside the sec- 
ond layer (S2). P2{si,S2) will then reduce to a product 
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Pi{si)Pi{s2), exactly as required for the multiplication 
rule to hold. Note that the Siegert relation (|l^ ) implies 
that the multiplication rule applies to g2{T) — 1 as well. 



V. EXPERIMENTAL SETUP 

Our experimental setup consists of a frequency- 
doubled NdYAG laser ("Verdi" from Coherent, wave- 
length Ao = 532 nm) with a beam width of roughly 1 
mm. Only minor differences are found upon expanding 
the beam to 7 mm in width (see also Ref. ) . The laser 
radiation is directed at the surface of the sample consist- 
ing of two cells, and the multiple-scattered, transmitted 
light is collected using a single mode optical fiber. The 
intensity of collected light is analyzed by a digital corre- 
lator (see Fig. ||). Apart from the complex structure of 
the sample, the described experimental setup represents a 
classical one for diffusing- wave spectroscopy experiments 
in transmission geometry |^-^ . We pay special attention 
to the preparation of the sample, which is a key feature of 
our experiments. The first cell (thickness Li) is filled by 
a colloidal gel, prepared from a destabilized solution of 
polystyrene spheres (diameter 170 nm) at a volume frac- 
tion up to 20% in a buoyancy-matching mixture of water 
and heavy water [|ll|,0. In our experiments, the gel 
serves as a model nonergodic medium where the average 
mean square displacement of the constituent particles is 
reasonably well described by the model of arrested subd- 
iffusive motion [see Eq. (H)]. In general, gel systems un- 
dergo a complicated temporal evolution (see, e.g., Refs. 



0,yHl5|j5 

paper 



I), which is not a subject of the present 



For our measurements, we wait after the sol-gel 
transition till the height of the plateau of g2{T) — 1 is 
of the order 0.5, and the gel properties remain constant 
over the measurement time of about 10-30 min. Wc note 
that experimentally 32 (t) — 1 does not exhibit a perfect 
plateau but rather a long stretching. This indicates a sec- 
ond, long-time decay of the autocorrelation function, not 
included in our simple theoretical model [Eq. (|3|)]. The 
characteristic time scale of this decay extends to minutes 
or hours [ [Tl] , |5^ , and hence is beyond the time window of 
our experiments. 

Two different realizations of the two-cell sample were 
used in the experiments reported below (samples A and 
B for brevity). While sample A is a model sample con- 
structed to test the validity of our approach, sample B is 
optimized to facilitate its experimental application (the 
multiplication rule holds for this sample). 

Sample A. The sample consists of three disc-shaped 
microscope cover plates (radius R = 7 mm, thickness 
A — 0.15 mm) and two ring-shaped teflon spacers. The 
first cell contains a colloidal gel (volume fraction <!> — 
7%, photon transport mean free path II = 53 /im). The 
second cell is filled with a suspension of Ti02 powder 
(Aldrich, particle diameter < 5 /xm) dispersed in pure 
glycerol (refractive index n = 1.47). 



Laser 

532nm 



Fiber 



Correlator 
+ PC 



PM2 



PMl 



FIG. 3. Experimental setup. A laser beam (Aq — 532 nm) 
is incident on a sandwich of two turbid layers. The light trans- 
mitted through the sandwich is detected with a mono-mode 
fiber, and is subsequently analyzed with a photomultiplier 
(PM) /digital correlator unit. 

The transport mean free path for this suspension is 
I2 = 41 /im. The single scattering correlation time 
corresponding to the suspension in the second cell is 
To = 23 ± 1 s, as determined by standard DWS mea- 
surements in backscattering geometry The thick- 
nesses Li and L2 of the layers were determined from 
DWS measurements in transmission geometry by filling 
them with a suspension of known I* and fitting the re- 
sulting intensity autocorrelation functions with Eq. (p^). 
We find Li = 1.75 mm and L2 = 1-9 mm. The opti- 
cal thicknesses of both layers are comparable and high 
enough (Li/ll ~ 33 L2/I2 ~ 46 > 1), and therefore 
we call this sample "symmetric" . A typical photon path 
length s in the sample is s L'^/l* « 4 x 10^ mm (with 
L ~ Li -I- L2 ~ 4 mm, I* = I* 40 //m), which is 
two orders of magnitude larger than the sample thickness. 
Hence, the scattering of light in our sample is essentially 
multiple, and a typical photon path is expected to cross 
the buried inter-layer interface many times. Estimation 
[ |72[ of the attenuation factor / gives / ~ 0.98 ~ 1, and 
hence the multiplication rule will only hold for the sam- 
ple A in a trivial situation of aiLi <^ 1 and q;2-/j2 ^ 1, 
as Li/ll ^ ^2/^2- 

Sample B. The sample consists of two equally thick 
light-scattering cells (Hellma, Li = L2 = 1 mm). The 
thickness of the glass wall between the turbid media in- 
side the cells is A ~ 2 mm. Such a thick wall leads 
to a significant photon leakage in transverse directions, 
and we put an aperture of radius R = 2.5 mm in between 
the two cells to prevent the photons leaked out of the cell 
from reaching the detector. This leads to some additional 
decrease of the total transmitted intensity, while ensures 
the multiplication rule, as discussed in Sec. IV. Indeed, 
for / ~ 0.46 estimated theoretically ||7^ the multiplica- 
tion rule shoidd hold whenever Ln/l^ S> 1/(1 — /^) ^ 1 
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and an^n ^ 1 — 1, which coincides with the condi- 
tions of vahdity of Eq. (|l^). We expect the actual value 
of / to be even smaller due to (total) reflections of light 
at the interface which are not included in our present 
theoretical model. Even though the conditions of valid- 
ity of the multiplication rule are satisfied, to apply it in 
a real experiment [i.e. to find g2^\Li,T) — 1 as a ratio 
of g^\Li, L2,t) — 1 and g2^\L2,T) — 1] we need both 

g^^^Li, L2,t) — 1 and g':^\L2,T) — 1 to be essentially 
different from 0, as otherwise experimental errors in de- 
termination of the above autocorrelation functions can be 
significant, making the application of the multiplication 
rule impractical. To overcome this problem, we choose 
the second light-scattering cell to be of much smaller op- 
tical thickness than the first one, while the motion of 
scatterers is made slow inside it. Namely, the first cell is 
filled with a very turbid gel (volume fraction $ = 20%, 
l\ — 19 /im), while the second one contains a moderately 
turbid suspension of colloidal polystyrene dispersed in 
a water-glycerol mixture {I2 = 300-500 /im). The de- 
crease of the measured autocorrelation function due to 
the motion of scatterers in the second, ergodic cell is 
thereby shifted to long times r ~ r2 = tq{12/L2Y and 
does not obscure the light-scattering signal of the first, 
nonergodic cell. Hence, the multiplication rule can be 
efficiently applied for t < T2, while for t ^ r2 the mea- 
sured autocorrelation function will decrease to due to 
the motion of scatterers in the second cell. For the sam- 
ple B, Li/l\ w 50 > ^2/^2 ~ 2-3, and therefore we call 
this sample "asymmetric" . 

VI. RESULTS AND DISCUSSION 

In order to demonstrate the efficiency of the method 
proposed and theoretically justified in Sees. and 
to show the feasibility of diffusing-wave spectroscopy in 
nonergodic media, we have carried out several model ex- 
periments using samples A (symmetric) and B (asym- 
metric) . 

A. Symmetric sample 

Open circles in Fig. ^(a) show the intensity autocor- 
relation function g\ {Li, L2,t) ~ 1 of light transmitted 
through the symmetric sample A. The solid line going 
through the open circles of Fig. ^(a) is obtained by fit- 
ting the experimental data with Eq. jll] ) assuming / = 1 
(no leakage of light at the wall separating the cells). The 
fitting parameters are the gel parameters S = 2.24 nm, 
Tc = 25 fis, p — 0.6 [see Eq. (|)], while the parame- 
ters of the medium in the second cell, as well as the cell 
thicknesses Li and L2 were determined independently. 
We note that the two sets of parameters are only weakly 
correlated in the fit. Theory and experimental data are 
found in excellent agreement. 




FIG. 4. (a) Intensity autocorrelation function for transmis- 
sion through the symmetric two-cell sample A (open circles), 
and through the isolated second cell (full squares') . Theoret- 
ical fits [Eq. (|l^ with / = 1 and Eq. (|l|)] are shown by 
solid lines. The dotted line shows the theoretical correlation 
function of the two-cell sample with the first layer assumed 
to be rigid (tc —> 00). (b) Different contributions to the in- 
tensity autocorrelation function are shown. The solid line is 
the fit to the two-cell data of the panel (a). To obtain the 
dashed-dotted line, we keep the second layer rigid (ro 00), 
while all other parameters are unchanged. The dotted line is 
that of the panel (a) multiplied by the plateau value 0.422. 

As shown in Fig. 0(a), the autocorrelation function of 

the two-cell setup 52 (^i: ^2, t) — 1 exhibits two charac- 
teristic decay times. The fast decay is due to the gel in 
the first layer, while the slow decay corresponds to the 
dynamics in the second layer. This is illustrated in Fig. 
^(b), where we decompose g^\Li, L2, r) — 1 in two con- 
tributions due to the dynamics in the first and second 
layers, respectively. The contributions of the first (dash- 
dotted line) and second (dotted line) layers are obtained 
by assuming the second or the first layer to be rigid [Eq. 
( [ll| ) with f — 1 and 02 = or ai = 0, respectively]. 
The contribution of the second layer is renormalized so 
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that its value at r = is equal to the contribution of the 
first layer at r — > cx). The figure suggests that a simple 
interpretation of the correlation function resulting from 
the two-cell setup in terms of contributions of individual 
layers can be given if the decay times due to different 
layers are well separated. 

The dotted line in Fig. 0(a) shows the contribution of 

(2) 

the second layer to {Li,L2, t) — 1 [the same as in Fig. 
^(b) but renormalized to 1] . Due to the presence of the 
first layer, although assumed to be rigid, the decay of 
the autocorrelation function becomes significantly faster 
as compared to the second layer taken alone [full squares 
in Fig. ^(a)]. This stems from the fact that the typical 
lengths of photon paths in the second layer are increased 
due to the presence of the first layer, which acts as an 
effective "diffuse mirror" , increasing the probability for 
a light path to get back to the second layer instead of 
being diffusely reflected from the sample. 

The measurements performed on the symmetric sam- 
ple A allows us to conclude that the multiple scattering 
of light in a two-cell sample is correctly described by the 
theory developed in Sec. [II. At the same time, the re- 
sults of this subsection serve as a test of our method, al- 
lowing diffusing-wave spectroscopy to be applied to non- 
ergodic random media. Indeed, the parameters of the 
(nonergodic) gel {5, Tc and p) can be obtained from the 
fit to the experimental data [see Fig. ^ , provided that the 
parameters of the medium in the second cell as well as 
the thicknesses of both cells are measured independently. 



B. Asymmetric sample 




Although the results of the previous subsection seem to 
be sufficiently convincing to justify our method of per- 
forming DWS in nonergodic media, we will now show 
that the interpretation of experimental data can be fur- 
ther simplified by introducing leakage (or absorption) of 
light at the interface between the layers in combination 
with a proper choice of parameters of the second layer. 
An optimal choice of the parameters corresponds to our 
asymmetric sample B (see Sec. For this sample, the 
optical thickness of the second layer is much smaller than 
that of the first one, and the leakage of light at the inter- 
layer wall is considerable. This ensures the validity of the 
multiplication rule formulated in Sec. as confirmed 
by the experimental results presented below. Figures 
H(a) and (b) show the intensity autocorrelation functions 

(2) 

52 {Li, L2,t) — 1 obtained for the two-cell setup (open 
circles) . The results for two different second cells are dis- 
played. Full squares show g^\L2, t) — 1 obtained for the 
isolated second cells, just as in Fig. ^. The correspond- 
ing single scattering correlation time tq = {k^Ds)'^ of 
the particles in the second cell has been varied by chang- 
ing the glycerol content in the suspending liquid, which 
affects the viscosity of the liquid and, consequently, the 
particle diffusion coefficient Db- 



1E-8 1E-6 1E-4 0.01 1 
X[S] 

FIG. 5. Multiplication rule: (a) Intensity autocorrelation 
function measured in transmission through the asymmetric 
two-cell sample B (open circles). Full squares show the auto- 
correlation function measured for the isolated second cell [de- 
cay time T2 — To(i2/i2)^]. (b) Same as (a), but for a different 
medium inside the second cell (higher glycerol content of the 
suspending liquid leading to a reduced particle diffusion coef- 
ficient and a larger value of r2). (c) The ratios of the autocor- 
relation functions are shown for the data of panels (a) (open 
circles) and (b) (full squares). The results are identical for 
T < 5t2 [with T2 of panel (a)], and yield the ensemble-averaged 
intensity autocorrelation function (?2^''(Li,r) — 1 correspond- 
ing to the nonergodic first cell taken alone. Inset: Relative 
difference A(r) = {[gi'^\r) - 1] - [gi'^r) - l]}/[ff^"'(r) - 1] 
between the two correlation functions of the main plot. 

According to the multiplication rule [Eq. (^6|)] , the ratio 
[g^\Li, L2, T) — l]/[g2^\L2, t) — 1] is expected to be equal 
to g2^\Li,T) — 1, i.e. the autocorrelation function cor- 
responding to the (first) nonergodic light-scattering cell. 
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To check this prediction, we have calculated correspond- 
ing ratios for the curves of Figs. 1(a) and (b). The results 
are presented in Fig. ^(c) [full squares correspond to the 
data of Fig. H(a), while open circles — to the data of Fig. 
1(b)] . The two curves are indistinguishable, supporting 
the validity of the multiplication rule. Experimentally, 
we find that for r < 5t2 [with T2 = tq{12/L2Y] the de- 
viations from the multiplication rule are negligible. The 
perfect agreement of the results obtained for two differ- 
ent media in the second cell confirms the validity of our 
method. 

As an additional check of the multiplication rule, we 
have performed measurements of the intensity autocor- 
relation function f/2 (^ii ^2, ''') — 1 with the first layer of 
the double-layer sample being totally rigid (white paper 
was used to model rigid but turbid random medium). 
The results are presented in Fig. |^ by open circles for 
f/j^-* (Li, L2, ''■) ~ 1 and full squares for (/2^-'(L2, t) — 1 (the 
latter is measured for the isolated second layer) . Full cir- 
cles are obtained by calculating the ratio [g^ {Li, L2,t) — 
M/[92^\-^2,t) — 1]. As expected, the ratio is one over 
a sufficiently extended time range (the ratio one corre- 
sponds to no scatterer motion in the paper). 



out some pitfalls when designing a two-cell experiment. 
Firstly, it is worthwhile to mention the problem of low 
transmitted light intensity. Building a sandwich with reg- 
ular glass cuvettes can result in a sample of significant 
overall thickness. Light incident on the setup spreads 
out diffusely which leads to relatively low intensity of 
light reaching the transmission side. The problem be- 
comes even more severe if absorption and/or leakage of 
light are introduced at the inter-cell wall to simplify the 
interpretation of experimental data by using the multi- 
plication rule. Secondly, due to reflections and additional 
leakage of diffuse light at the interface glass walls, some 
of the light scattered out of the cell can reach the detector 
without passing through the whole sample. To avoid this 
problem, we recommend to put an aperture (e.g. a black 
plastic foil with a centered hole of typically 5 mm in di- 
ameter) between the cells to suppress photons scattered 
snake-like along the cell walls. Ideally, the double-cell 
sandwich would have only one moderately absorbing thin 
interface wall with a built-in circular aperture forcing all 
transmitted light to pass through the whole sample. 

VIII. SUMMARY AND OUTLOOK 




1E-7 1E-6 1E-5 1E-4 1E-3 0.01 0.1 1 

t[s] 

FIG. 6. Intensity autocorrelation function for a two-cell 
sample with a rigid first layer (the colloidal gel in the first 
layer was replaced by a white paper which mimics a rigid tur- 
bid medium) (open circles). Full squares show the autocor- 
relation function corresponding to the isolated second layer. 
The ratio of the two autocorrelation functions (full circles) is 
one to a good accuracy for t < T2 = to (/1 7^2)^, correspond- 
ing to the absence of scatterer motion in the first layer. 



VII. HOW TO APPLY THE TWO-CELL 
TECHNIQUE ? 

One of the major advantages of the two-cell technique 
is its simplicity. Nevertheless, we would like to point 



In the present paper we propose and test experimen- 
tally a new method for the application of diffusing-wave 
spectroscopy to nonergodic turbid media. We show that 
light transmitted trough a sandwich of two turbid sam- 
ples can be considered ergodic even if only the second 
sample is ergodic. The autocorrelation function of the 
transmitted intensity can be quantitatively described by 
the diffusion theory. This allows direct application of 
diffusing-wave spectroscopy for the characterization of 
nonergodic media without any additional efforts usu- 
ally required to achieve ensemble averaging (i.e. with- 
out translation or rotation of the sample in course of the 
correlation function measurement [ p9| , p4 35 5^]). Proper 
averaging of the light-scattering signal is ensured by 
adding a second, ergodic light-scattering cell with un- 
precedented accuracy. In order to simplify the analysis 
of the experimental data, the parameters of the double- 
cell sample can be optimized. Namely, moderate absorp- 
tion and/or leakage of light should be introduced at the 
interface between the light-scattering cells, and the op- 
tical thickness of the second, ergodic cell should be re- 
duced well below the optical thickness of the cell con- 
taining the nonergodic medium, while the dynamics of 
scatterers in the second cell should be chosen slow. Un- 
der these conditions we have shown that the field au- 
tocorrelation function of light transmitted through the 
double-layer sample can be written as a product of au- 
tocorrelation functions corresponding to the individual 
layers: g^j^\Li, L2,t) ^ g[^\Li,T) x g[^\L2,T). Conse- 
quently, it is sufficient to measure the intensity autocor- 
relation functions [^^ of the two cell setup g^\Li, L2, t) 
and of the second, ergodic layer gj^'' (i2, '''), in order 
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to obtain the properly averaged intensity autocorrela- 
tion function of the nonergodic layer (72^''(Li,r) — 1 = 



[(?2^''(Li,L2,t) — \\/[g'^' {L2,t) — 1]. Henceforth applica- 
tion of diffusing-wave spectroscopy for the characteriza- 
tion of nonergodic media becomes straightforward. Our 
method is particularly suited for the study of time evolv- 
ing systems, e.g. aggregating and gelling particle sus- 
pensions (see also Refs. Q and [^), since the data ac- 
quisition time can easily be adjusted from a few seconds 
to several hours. Other promising applications include 
DWS-microrheology which has previously been limited 
to time and length scales where the motion of tracer par- 
ticles is not significantly constrained |^,^. We expect 
the two-cell technique to extend the measurement range 
of this method and to provide access to a novel class of 
(solid-like) materials. 
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